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Abstract 

Using methods of effective field theory, we derive the first all-order factorization theorem 
for the Higgs-boson production cross section with a jet veto, imposed by means of a stan- 
dard sequential recombination jet algorithm. Like in the case of small-g^ resummation in 
Drell-Yan and Higgs production, the factorization is affected by a collinear anomaly. Our 
analysis provides the basis for a systematic resummation of large logarithms ln(m# /p™ to ) 
beyond leading-logarithmic order. Specifically, we present predictions for the resummed 
jet-veto cross section and efficiency at next-to-next-to-leading logarithmic order. Our 
results have important implications for Higgs-boson searches at the LHC, where a jet 
veto is required to suppress background events. 



1 Introduction 



While the strong interactions generally complicate physics at hadron colliders, the Higgs- 
boson production cross section in gluon-gluon fusion is one of the few places where QCD is 
trying to help: radiative corrections substantially increase the production rate at both next- 
to-leading (NLO) [U |2] and next-to-next-to leading order (NNLO) in perturbation theory 
[3j IH [5] . However, to interpret results from Higgs searches and to study Higgs properties, it is 
important that the higher-order corrections are under control. In [6j [7] we have analyzed the 
total production cross section and argued that the largest effects arise in the virtual corrections 
to the gg — >■ H amplitude, which are given by the scalar gluon form factor evaluated at time- 
like momentum transfer. These corrections can be resummed using renormalization-group 
(RG) methods. 

The total cross section is not directly relevant for Higgs-boson searches, since these impose 
a jet veto fij, < p™ to to enhance the signal-to-background ratio. In particular, the jet veto is 
essential to reduce the background to the pp — > H — > W + W~ process arising from top-quark 
production with subsequent t — > bl + u decay. To get sufficient background reduction, the 
experimental searches reject all events containing jets with transverse momenta larger than 
p™ to 20-30 GeV. It was observed that the i^-factor for Higgs production gets reduced when 
such a cut is imposed [H [91 [10]. Keeping the parton distribution functions (PDFs) and the 
value of a s fixed, the total cross section increases by more than 100% when going from LO to 
NLO, and by another 35% when going from NLO to NNLO. The corrections to the jet-veto 
cross section are however much smaller. This is evident from the first two plots in Figure [U 
which show fixed-order results for the jet-veto cross section for Higgs-boson production at the 
LHC, obtained using different scale choices. We employ a standard jet-clustering algorithm 
with R = 0.4 (see below) and use MSTW2008NNLO PDFs [H]. For example, at $5 to = 20 GeV 
the difference between the NLO and NNLO cross sections is only about 10%. Note that in 
the presence of the jet veto the cross section involves the two disparate scales, p^ to and m#, 
and it is not clear a priori which value of the factorization and renormalization scales one 
should choose (we set p = pj = p r throughout). Indeed, the results look rather different 
depending on whether one takes p = m H (first plot) or p = p^ ta (second plot) as the default 
value. It thus appears somewhat suspicious that the NNLO results for p™ to > 15 GeV happen 
to be rather similar in the two cases. In [7] we have argued that the reduction of the i^-factor 
observed in the presence of a jet veto arises from an accidental cancellation of the large positive 
corrections to the total cross section with large negative corrections from Sudakov logarithms 
a™ ln 2n (m#/p™ to ), which arise because the jet veto enhances the contributions from emissions 
collinear to the beams. 

The presence of large Sudakov logarithms becomes obvious when one considers the jet- veto 
efficiency e(p™ to ) = a(p^ to )/a to t, defined as the ratio of the vetoed to the total cross section. 
In this quantity the common large virtual corrections to the gg — > H amplitude drop out, 
leaving only the large Sudakov logarithms behind. The third panel in Figure [1] shows that 
the veto efficiency indeed receives very large higher-order corrections. Even at NNLO the 
remaining scale variation is rather significant. In the figure we choose a high default value for 
the factorization scale. Choosing p ~ p™ to would lead to a complete breakdown of the fixed- 
order expansion for the efficiency. In addition to the large scale uncertainty the result suffers 
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Figure 1: The jet-veto cross section a(p™ to ) and efficiency e(p™ to ) for Higgs production at 
the LHC, at NLO (shaded light bands) and NNLO (dark bands) in fixed-order perturbation 
theory. The first two panels compare the scale choices fi ~ m# and /i ~ p^ to . The right 
panels show predictions for the veto efficiency, defined either as a perturbative series in a s 
(third plot), or as the ratio of the n th -order vetoed cross section to the n th -order total cross 
section (last plot). 



from a significant scheme ambiguity. The last plot in the figure shows the predictions for the 
efficiency obtained when one expands the two cross sections a(p™ to ) and <T to t to NLO or NNLO, 
but then takes their ratio without performing a further expansion in powers of a s . Comparing 
the last two panels, we observe that the scheme dependence remains uncomfortably large even 
at NNLO. 

The leading-logarithmic (LL) corrections to the jet-veto cross section were studied us- 
ing parton showers [T2"| [13] , but for a long time no systematic resummation of higher-order 
logarithmic terms was available. To improve the accuracy of the numerical predictions, the 
parton-shower results were re-weighted to the Higgs-boson spectrum obtained at next- 
to-next-to- leading logarithmic (NNLL) order, as implemented in the code HqT [TJ]. As an 
alternative, it was suggested to use the event-shape variable beam thrust instead of a jet veto 
to discriminate against the background from top decays [15]. The NNLL resummed results for 
this variable confirmed the picture that there are large corrections from collinear emissions. 
It was found that a fixed-order computation is not reliable and that the scale variation un- 
derestimates the perturbative uncertainties. While beam thrust is theoretically simpler than 
the jet veto, it would be more difficult to use in experimental analyses, since this observable 
is sensitive to underlying-event and pile- up effects. Based on the results for beam thrust, 
an alternative way to estimate the scale uncertainties for Higgs searches using jet bins was 
proposed in [To] . 

Very recently, however, it was pointed out that the resummation for the jet-veto cross 
section defined with a standard sequential jet algorithm is indeed feasible, and numerical 
results for the cross section at NLL order were presented [T7j. These results were obtained 
using the Caesar code |18j . which performs automated NLL resummations for observables 
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that are subject to certain conditions, such as recursive infrared safety and globalness. In 
this paper, we study the resummation for the cross section in the presence of a jet veto in 
Soft-Collinear Effective Theory (SCET) [HJ |20j EJ [22]. In Section El we propose an all- 
order factorization theorem, according to which the cross section factorizes into a product 
of a hard function, a soft function, and two beam-jet functions. This simple factorization 
into a product of functions (rather than a convolution) arises because, at leading order in an 
expansion in powers of p^ to /mH, the &T-type jet algorithms do not cluster soft and collinear 
radiation inside the same jet. As for small-gr resummation, the naive factorization is affected 
by a collinear anomaly [23], which induces a dependence on the Higgs-boson mass in the 
product of beam-jet and soft functions. We derive the all-order form of this anomaly and 
present simple analytic formulae for the resummed cross section. In Section [31 we then give all 
ingredients required for carrying out the resummation at NNLL order, including a certain two- 
loop coefficient capturing the dependence of the resummed cross section on the jet algorithm. 
We extract this coefficient from partial results at NNLL order presented in the appendix of 
|17j . Phenomenological predictions for the jet- veto cross section and efficiency are presented 
in Section [H In the final section we summarize our findings. 



2 Factorization and resummation of the cross section 
2.1 Preliminaries 

In the heavy top-quark limit, the effective Lagrangian describing Higgs production via gluon- 
gluon fusion reads 



^^k^-g;^, (i) 

1Z7T V 



where the Wilson coefficient Ct accounts for higher-order loop effects. At NLO in RG-improved 
perturbation theory, it is given by |7] 

Ct{m t ,(i) = l+ ^ (5C A - 3C F ) + — — + . . . , (2) 

where fi t ~ irk is the matching scale at which the top quark is integrated out, and /?q = 11— § ft/ 
and j3i = 102 — y rif (for N c = 3 colors) are the first two expansion coefficients of the QCD 
/3-function. In the presence of a jet veto, the differential cross section for Higgs production at 
the LHC can then be written as 

da( P v f°) = — ( ^MVc2 (m 2 ) <Pq i d ± iq . x 

yFT ' 2s\12irvJ tK 1 ^ 2vr 3 2E„ J 

V / (3) 

x (P( Pl )P(p 2 )\G; v G^(x) \X) (X\ G b pa G^ b (0) \P( Pl )P( P2 )) , 
x 

where the prime on the sum indicates that we only sum over those hadronic final states X that 
pass the jet-veto cut. We work with the usual class of sequential recombination jet algorithms 
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|25j . with distance measure 



dij = imiL{pTi,PTj) 



R 



(4) 



diB — PTi 



where n = 1 corresponds to the fcy algorithm [2S1 EZ], n = to the Cambridge- Aachen 
algorithm [281 [29], an d n = — 1 to the anti-fc-r algorithm [30]. The particles with the smallest 
distance are combined into a new "particle", whose momentum is the sum of the momenta 
of the parent particles. If the smallest distance is diB, then particle i is considered a jet and 
removed from the list. The procedure is iterated until all particles are grouped into jets, 
i.e., the algorithm is inclusive. The corrections to the cross section which are logarithmically 
enhanced at small p^ to arise from emissions that are soft or collinear to the proton beams. 
In perturbation theory, a sensitivity to the jet algorithm appears first at NNLO through 
diagrams with two real emissions. At this order a dependence on R arises, but the cross 
section is still independent of n. The two particles are clustered into a single jet whenever 
a/ Ay 2 2 + A0 2 2 < R, and otherwise they are treated as two jets. 

We will analyze the cross section using the formalism of SCET, in which highly energetic 
particles aligned with the colliding protons are described in terms of collinear and anti-collinear 
quark and gluon fields, and soft particles emitted from the beam jets are described in terms of 
soft fields. The effective theory implements an expansion of scattering amplitudes in powers 
of the small parameter A ~ p v rf to /mH, where the jet veto sets the characteristic size of all 
transverse momenta in the process. We introduce two light-like reference vectors n M and W 1 
(satisfying n ■ n — 2) parallel to the beam axis and decompose all 4-vectors in the light-cone 
basis spanned by these vectors, 



The different types of modes relevant to our discussion are characterized by the scalings of 
their momenta (p + ,p_,p±) with powers of A, namely p^ ~ m#(A 2 , 1, A) for collinear particles, 
Pg ~ m#(l, A 2 , A) for anti-collinear particles, and p^ ~ mjy(A, A, A) for soft particles. Hence, 
the particles in these three categories have transverse momenta of order the jet veto, but very 
different rapidities. 

At leading order in power counting, the gluon operator in ([I]) can be matched onto an 
operator in SCET consisting of gauge-invariant, collinear and anti-collinear gluon fields A£ 
and *4g. The formalism of expressing SCET operators in terms of gauge- invariant building 
blocks is described, e.g., in [2H HD]. In the present case, the corresponding matching relation 
reads [7] 



where q 2 = m 2 H is the time-like momentum transfer carried by the current. The explicit 
expression for the hard matching coefficient Cs at NLO in RG-improved perturbation theory 
is given in relation (IA4[) of the appendix. Note that only the physical (transverse) polarization 
states of the gluons appear in the above relation. The objects 5* n and are soft Wilson lines 



p^ = n ■ p — + n ■ p — + p^_ = p+ + pt + p^_ ■ 



(5) 
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in the directions of the light-like vectors n and n, which extend from the interaction point to 
infinity. Assuming that all quantum fields vanish at infinity, the product S^Sn corresponds 
to a gauge- invariant Wilson loop with a cusp at the point where the two Wilson lines come 
together. In ([6]) we have already applied the SCET decoupling transformation [20], which 
removes the interactions between (anti-)collinear and soft fields from the SCET Lagrangian. 
It is then always possible to group the fields of different types together and factorize any SCET 
operator into gauge-invariant products of collinear, anti-collinear, and soft fields. Likewise, 
the final states \X) can be written in the factorized form \X) = \X C ) cg> \X 5 ) (g) \X S ). 

2.2 Jet clustering and mode separation 

We now exploit the fact that the jet clustering algorithm does not group particles with different 
momentum scalings (collinear, anti-collinear, or soft) into the same jet. The reason is that, 
generically, the rapidity difference between two such particles is of order ln(m#/p™ to ), which 
is parametrically larger than the radius parameter R employed in the jet algorithm (J4]). Only 
in corners of the phase space, e.g. when a soft emission becomes collinear to the beam, soft 
and collinear radiation can be combined. However, since the cross section does not exhibit 
additional singularities in the corresponding limit, such configurations only give rise to power- 
suppressed contributions. 

To make this argument more precise, we must first specify the power counting adapted for 
the radius parameter R, which we choose such that 

A < R < In - . (7) 
A 

If R was too large, the argument about clustering given above (and explained in more detail 
below) would fail, while for too small R large logarithms ln n i? would arise, which would 
require a special treatment. These "clustering logarithms" have a complicated structure in 
higher orders [521133] , and it is currently not understood how to resum them. Assuming that 
R satisfies relation ([7]), we will derive a factorization theorem valid at leading power in A in 
the formal limit where A 1. Once the result has been established, it will hold for any value 
R = 0(1). Even though A is not very small in practice, the above condition is satisfied for 
realistic values R = 0.4 or 0.5 used in experimental analyses. For instance, with A = p^ to /rriH 
we have A ~ 0.16 and ln(l/A) ~ 1.83 for p™ to = 20GeV and m# = 125 GeV. Our resummed 
results for the cross section will receive power corrections in A, which will be added when we 
match them to fixed-order perturbation theory. 

Consider now the behavior of soft, collinear, and anti-collinear modes under the jet clus- 
tering algorithm (jl]). For two modes from the same sector, the quantities and diB are 
parametrically of the same order, and depending on the precise values of the transverse mo- 
menta, rapidities, and polar angles of the two particles, they are either clustered together or 
considered two jets. On the other hand, for two modes from different sectors (e.g. a soft and 
a collinear particle) we have parametrically dij ~ A n ln(l/A) but diB ~ A n , where it is crucial 
that R = 0(1) or smaller. It follows that generically the two particles are not clustered into 
the same jet. Since the soft and (anti-) collinear modes have the same virtuality, they live 
along a hyperbola in the (p + ,p_) plane, and their precise separation along this hyperbola is to 
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some extent arbitrary. The fact that these modes differ by large rapidities gives rise to large 
logarithms, which are accounted for by the collinear anomaly. In complete analogy with the 
construction of the SCET Lagrangian, where based on the generic scalings of the fields one 
does not include soft-collinear interaction terms, it is unnecessary to consider the degenerate 
case where a collinear and a soft mode near the boundary are clustered into a single jet. Since 
there are no enhancements of the cross section in these power-suppressed phase-space regions, 
boundary effects do not contribute at leading power. 

One may also worry about the contributions from modes with smaller virtualities, p 2 <C 
(p™ to ) 2 . For example, an on-shell soft mode, which accidentally is closely aligned with the 
beam axis, would have momentum scaling (A 2 , A, A 3 / 2 ). This mode has the same rapidity as 
a collinear mode, and hence it may be clustered with such a mode. Indeed, it may also be 
regarded as a collinear mode, whose minus component is accidentally small. The important 
point is that, because of their small transverse momenta, such modes play no role for the total 
Pt of a jet. Therefore, an arbitrary number of them can be emitted, and their effect cancels 
out in the factorization theorem. This is in analogy with the cancellation of ultra-soft modes in 
the factorization theorem for the Drell-Yan cross section at small transverse momentum [23] . 

It follows from this discussion that, at leading power in A, the cross section factorizes into 
the production of soft and (anti-)collinear jets, and the jet veto can be applied separately in 
each sector. Changing variables from d 3 q to dyd 2 q±, where q± and y denote the transverse 
momentum and rapidity of the Higgs boson, we obtain from fl3]) 

d °(PT°) = ±(^) 2 CUnil»)in,' H \C s (- m %, l i)\> ggi J A«* 

(P(pi)| Ar(x + + x ± ) \X C ) {X c \ A»/(0) |P(pO> g L „gi 



X 



l 1 



X 
X 



(P(p 2 )\A p /(x. + x ± ) \X B ) (X B \ Ar(0) \P(p 2 )) 

Xc 

( I (SiS n ) a \x ± ) \X S ) (X s \ {SlS n ) cd (0) | ) + . . . , 



(8) 



X 

Xc 



where the dots represent power-suppressed terms in the SCET expansion parameter A ~ 
p™ to /mH, which we neglect. Above, we have already performed the multipole expansion of 
the SCET fields, which neglects those components of x in the arguments of the fields that 
would introduce subleading momentum components [22J. To identify these components we 
have used that ~ (1, 1, A -1 ), which is conjugate to the momentum g M ~ m#(l, 1, A) of the 
Higgs boson. 

We now use color conservation for the individual matrix elements and define 

B^(z,x ± ,pT,») = f dte-^^(P(p)\A^(tn + x ± ) \X C ) (X c \ ^(0) \P{p)) , 

x c 

S(x ± ,?r^) = (0 I (^) a6 (* ± ) I**) (X.\ (SlS n ) ba (0) |0> . (9) 
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(10) 



The function B% v describes the structure of the jet of collinear particles inside one of the 
colliding protons (the one moving along the light-like direction n M ), which is probed at small 
transverse distances x±. In the context of SCET, such functions are referred to as beam 
functions [M]. The corresponding function B^ u for the other beam jet is given by the same 
formula with the replacements n — > n and c — > c. The soft function S describes the physics 
of soft gluons emitted from the colliding beam particles. Note that both the beam functions 
and the soft function depend on the jet algorithm used for the jet veto. This is implicit in our 
notation, which only makes the dependence on p™ to explicit. 

When inverting the first relation in ([9]), one must take into account that the beam function 
B£ u (z, x±,p^ to , /x) has support for — 1 < z < 1, and it can be shown that B£ u (—z, x±,p v rf to , /S) = 
—B^ u (z, — xj_,p v rf to , /i) (see e.g. (35]). Restricting the integration range to positive z values thus 
leads to an extra factor of 2 for each beam function. After some straightforward algebra, we 
obtain 

da{pT) = vMC 2 t {m 2 t ^) \C s (-m 2 H ,fi)\ 2 ™* dy p± [ d 2 x ± e~^ 

x 22^(6, *±,P V T°, v) x ± ,^ cto , //) <S(x ± ,/ T eto , aO , 

where £1,2 = \fr ' e ±y and r = {m 2 H + The Born-level cross section is 

m 2 H aftp) 

Up to corrections of order A 2 , one can replace r — > m 2 H /s in (flQl) . Given that the peak of the 
transverse- momentum distribution is around 10 GeV [36], this is a very good approximation. 
Integrating over q±, we then obtain the simple result 

dCX ^ T - = &o(y)C?( m t,y) \Cs(-m 2 H ,fi)\ 2 

dy (12) 

x w? (6, o, pT, aO o, pT, h) 5(0, pt°i ») ■ 

Since the jet veto does not prefer any direction in the transverse plane, we can replace 

£f M,/ T ct >) = ^-B c (z,pT°,y), (13) 



where 



B c (z,pT°,») = -^ [ ^e-^^^P(p)|^i a (tn)|X c >(X c |^ ±/i (0)|P(p)>. (14) 

Note the similarity of this function with the usual definition of the gluon distribution function 
in the proton, which in SCET notation reads 



^/p(*,A0 = -^ / dte-^(P(p)\A^(tn)A a c±fl (0)\P(p)}. (15) 



Dropping the vanishing argument x± = in the soft function, we finally obtain 

dy 



} =a (^C 2 (m 2 , f i) \C s {-m 2 H ^)\ 2 B c (Ci,pT°, y)B- c {t 2 ^ to ^)S{p^°, fi. (16) 
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2.3 Collinear factorization anomaly 

The above formula appears to disentangle the hard scales m t and m# from the soft scale 
pveto ^ rriH- The factorization is however not yet complete. Closer inspection shows that the 
individual ingredients to the cross section (TIB]) are ill-defined in dimensional regularization, 
so that an additional regularization is needed in an intermediate step. The situation closely 
resembles the case of the Drell-Yan cross section at small transverse momentum analyzed in 
[23] . As in this case, it is convenient to employ an analytic regulator for this purpose. 

It was recently shown that the additional analytic regularization is only needed for the 
phase space, and can be achieved by replacing the usual phase-space integration by [37] 

Jd d k5(k 2 )6(k°) Jd d k^-^jY 5(k 2 )6(k°). (17) 

The regularization softens light-cone singularities in the cross section. It introduces a new 
scale v, which plays an analogous role to the scale /i entering in dimensional regularization. In 
the above form, the analytic regularization leaves the structure of the effective theory intact, 
and all steps in the derivation of the factorization formula ffTBT) remain the same, except that 
the phase-space integrals implicit in the sums over intermediate states Yl'x ^ n ® must be 
regularized according to (jT7|) . Note that the symmetry n -H- n between the collinear and 
anti-collinear beam-jet functions is broken in the presence of the regulator, which involves the 
n ■ k component in both sectors. An alternative analytic regularization scheme, based on a 
regularization of the SCET Wilson lines, was proposed in [3BJ EH]- 

In the presence of the extra regulator, the jet and soft functions are separately well de- 
fined, but individually they contain singularities (l/a n poles) in the analytic regulator. These 
cancel in the product of the jet and soft functions entering a physical cross section, so that the 
regulator can be removed at the end of the analysis. Importantly, however, the a regulator 
breaks the invariance of the soft and jet functions under a rescaling of the large light-cone 
components of the momenta of the external particles. This symmetry of the low-energy theory 
is anomalous: it is not recovered after the regulator is removed. Instead, one finds an anoma- 
lous dependence of the low-energy theory on the hard scale m#. Following j23j HQ], we now 
derive the all-order structure of this dependence by solving differential equations expressing 
the independence of the cross section of the analytic regularization scale v. 

In order to constrain the anomaly, we first slightly generalize our observable. We veto an 
event if it contains left-moving jets with transverse momentum above a value p v £ to or right- 
moving jets with transverse momentum above p v ^ to . The beam-jet functions are the same 
as before, with p^ to replaced by p v f to {p v r°) for the left-moving (right-moving) jet. The soft 
function, on the other hand, now depends both on p v £ to and p v ^ to , since soft jets can be emitted 
in any direction. Other than this, the form of the cross section is exactly the same as above. 
Our previous result ( FT6|) is recovered in the special case where p v £ to = p v ^ to = p^ to . Choosing 
the values differently will, however, allow us to derive additional constraints on the form of 
the anomaly. 

Since the original QCD expression does not require analytic regularization, we are 
guaranteed that the limit a — > can be taken in the product of the jet and soft functions in 
(I16j) . However, upon expanding in the regulator, the individual functions involve logarithms 



S 



of the scale v. Taking the logarithm of the product P = B C B 5 S of jet and soft functions in 
([TED, we have 



In P = In B c (in £1, pt°, + In B- c (in £ 2 ,pg°, /i) + In <S (in ^ pf , p™ to , /i) . (18) 

For /i ~ p™ to ~ PR to , no large logarithms beyond the ones we display explicitly can arise. Their 
functional form reflects the scaling of n ■ k in the various momentum regions, which was given 
after ([5]). The structure of this result closely resembles that encountered in our analysis of jet 
broadening presented in jlQ] . The fact that In P as well as its derivatives with respect to m# 
must be independent of v gives rise to a set of simple differential equations, which imply that 
the three component functions can be at most quadratic in their first argument, such that 

In^ln^^V) = / 2 (i) (6,^» ln 2 ^ 

^ ^ (19) 

+ /i ( %i,p V L Ct >) In^f- + / ( %i,p V L Ct » , 

and analogously for the right beam-jet function and the soft function, with expansion coeffi- 
cients and ft and logarithms ln(i//m#) and ln(z///i), respectively. The coefficients of 
the logarithmic terms are furthermore subject to the constraints 

f[ L) (6 , pT, /*) + (6 , $f °, n) + A (5) (Pl et °, pS to , /x) = o , 

1 (2°) 
/ 2 (L) (6,Pl cto ,/i) = fSPfafiP,!*) = - J fPifiFiPR ,!*) ■ 

Since the arguments in the various functions are different, it follows that f[ L ' RS) and f% must 
be independent of £i and £2- Moreover, the second constraint implies that the coefficients of 
the quadratic terms can only be a function of /i, 

fi L \ti,P?°,ri = f2(M)- (21) 
For the first-order coefficients, the symmetry of InP under L -H- R and £i -H- £2 implies that 

A (i?) (6,P v r^) - fi L \ei,p?°,ri = hpT^) + /i(pr» , (22) 

and 

& l \ZupT,v) + ti R \&pT^) + f^(pT°,p v R°^) 

(23) 

= fo(ti,pr o ,») + /o(6,^ eto ,M) + / (5) br to ,^ eto ^) • 

This last relation expresses that the dependence on £1 and £2 is governed by a single function 
/o, and a possible difference in the functional forms of and /g^ can be absorbed into the 

(S) 

soft function / . Using these results, we obtain 

InP = 2/2 (/i) In 2 ^ - In ™£ [^(p-*, /x) + /i(^ to , /*)] 

P (24) 

+ /o(ei,p V L Ct °,/i) + f fa,pT,H) + f^(p V L t0 ,PR t0 ^) ■ 
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Since the cross section is scale independent, the derivative of the above expression with 
respect to the renormalization scale must be equal and opposite to the anomalous dimension 
of the hard function, which is given in the appendix. This implies 

InP = -< US M In ^ - 4 7 %) , (25) 



din// cuspvrv \i 

where r^ sp is the cusp anomalous dimension in the adjoint representation, and 27 s = 7 5 + 
7* + j3(a s )/a s is the anomalous dimension of the gluon field as defined in [H]. Since this 
right-hand side of (J25]) is linear in ln(%/ju), we must require that 

/ 2 (A*) = 0. (26) 



But /2(a0 vanishes at tree level, and this equation then implies that it is zero to all orders 
in perturbation theory, as there is no way to compensate the \i dependence of the coupling 
constant. 

The outcome of the above considerations is that the logarithm of the product of beam-jet 
and soft functions is linear in ln(m#//x). This result is important for the resummation of 
logarithms. It implies that in the case p v £ to = p V R to = Pr to the product can be written in the 
form 



S c (ln^;ei,pr to ,/i) B- c (\n^-;&,pt ,») S (in ^ p? c P T > 



J q 2 =m? H 



VTiu \ I ma 



(27) 



^vcto / \ ^vcto 

Pl J \FR 



B g (^,pt , A*) B g fa, pf°, h) S g (pt°, PT, aO 



where we have defined 



F gg {pt°^) = h{pT^) 



ln£? g (£,0) = h(pt°^) In ^3 + /o(£,0) , (28) 

Pt 

In S g (pt , ft: H) = fo S) (PT, PT, /*) • 

The run dependence due to the collinear anomaly is now explicit in (1271) . Note also that the 
dependence on the longitudinal variables £1 and £2 is factorized and given in terms of a single 
function £> 5 (£,p™ to , p). The definition of these functions is not unique, since it is possible to 
shift p™ to -dependent terms into the soft function. For the special case p v f to = p v ^ to = p™ to , it 
is therefore natural to absorb the square root of the soft function into the definition of the 
beam-jet function B g . 

The structure of the result for the cross section, obtained when inserting (|2~T|) into is 
particularly simple for the case of the jet veto. In previous applications of the collinear anomaly 
in the context of small-gr resummation for Drell-Yan production [23] and jet broadening 
in e + e~ annihilations [ID], the anomalous exponents analogous to F gg depended on some 
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convolution variables, which were shared with the jet and soft functions in the factorization 
theorem. After integration over these variables, the resulting dependence on the hard scale 
can be quite complicated (see in particular the discussion in 02])- In the present case, on the 
other hand, the anomaly leads to an extra term in the cross section which, at fixed p™ to , is a 
pure power of m^. 



2.4 Short-distance expansion of the beam-jet function 

For the final factorization step, we now use that in practice the jet veto is much larger than the 
scale Aqcd governing non-perturbative hadronic effects in QCD. It is thus possible to calculate 
the physics associated with p^ to in perturbation theory, and to relate the refactorized beam-jet 



function B g S g to conventional PDFs [3H H3]. We write this relation in the form 



yl/2 

iju uuiivciiuiuiicii i Lji. -a yo^±v 1^0] ■ vvi 

ml dz 



f dz 

B g ^ to ^)s l J\ P t ,vt°^)= E / - W^>) <M£/^) 

„• z 



(29) 



i=g,q,q 



which is valid up to hadronic corrections suppressed by powers of Aqcd/p™* - By means of 
()27p and ( I29p . the cross section in ( II 6p can now be expressed in the final form 

- * (n)CKml»)\Cs(-m 2 H ,ri\ 2 (S) ( 3 °) 



dy w n t,r " Jy ""^ \,Pt 

' dz\ f 1 dz2 



r 1 dz i r dz2 



X 



where the sums over i, j extend over all flavors of partons (gluons, quarks, and anti-quarks). 
In the above factorization theorem, the dependence on the two disparate scales m# and 
p£ to <C rriH is completely explicit. Large logarithms of their ratio can be resummed by 
choosing the factorization scale around the value of the jet veto, /x ~ p™ to . With such a 
choice, the functions F gg and I g< _i have well-behaved perturbative expansions free of large 
logarithms. Explicit expressions for these functions will be presented in the next section. 
Large logarithms of m#/p™ to are contained in the anomaly term (in exponentiated form) and 
in the short- distance matching coefficients C t and C$, whose expressions are known at NNLO 
in RG-improved perturbation theory. They are collected in the appendix. 

A particularly nice formula is obtained after we integrate the differential cross section over 
rapidity. It reads 



°{pt°) = <r (n)Ct(mi,n) \C 8 (-m 2 B ,ti)\ 



^vcto 

ft 



-2F gg (p v f°^) 



(31 > 



f dz 

■ Jt Z 



X 

*>J=9)5 

where the Mellin convolution = I g <_t ® I g< _j of two kernel functions is defined as 

E ij {z,^°, l i) = / —I g ^(u,p^ : fi)I g ^(z/u : p^fi); i,j = g,q t (32) 

J z U 
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while ffij are flavor-singlet parton luminosity functions, given by ff gg = <p g /p ® 4> g /p, ff gq = 

To complete our discussion, we present the RG equations obeyed by the various component 
functions in (1301) . The evolution equation for the Wilson coefficient C t was first given in pM] . 
while that for Cs can be found, e.g., in [7]. For completeness, we list the corresponding RG 
equations in the appendix. The evolution equations for the functions F gg and I g ^ then follow 
from the scale invariance of the cross section. In analogy with our discussion in [23], we obtain 



F gg { P T^) = ^LM 



<iln/i 

^ v^p^v) = [2rl» in-^-2 7 %)] 

J - M 

where Pj^i are the usual DGLAP splitting functions. 



(33) 



3 Resummation at NNLL order 



For a consistent treatment at NLO in RG-improved perturbation theory, we require the 
anomaly exponent F gg at two-loop order, while the one-loop expressions for the kernel func- 
tions I g ^-i suffice. Up to 0{pt1), the general solution to the first evolution equation in (133]) 
reads 



ofio 



2 



d v 2 cto 



2 In 



^veto ' 
FT 



(34) 



where a s = a a (p)/(4m), Tq = 4Ca and T± = - *£ ) C\ - f C A T F n f are the one- and 
two- loop coefficients of the cusp anomalous dimension, and (3q = y Ca — § TpUf. The constant 
terms <i™ to and d™ to are independent of the jet veto p^ to but in general can depend on the jet 
algorithm. At one-loop order, the solution to the second evolution equation in ( 13"3"|) can be 
written as |23] 



I g ^ i (z,^ ,fj l ) = 6(l-z)6 gi 
where 7^ = -/3 , 

P«» = 8C A 

V { g l q {z) = AC F 



v- 

1 n 



7„ 9 £ 



7L g ^i(z) 



(35) 



1 - z 



L(i-*)+ 

l + (l-z) 2 



+ z(l-z) 



+ 2/? 5(l-z), 



(36) 



12 





mm® ®rrax? 




Figure 2: One-loop diagrams contributing to the anomaly exponent F gg and the matching 
coefficients I g ^ g (top row) and I g ^ q (bottom row). Vertical lines indicate cut propagators. 



are the one- loop DGLAP splitting functions, and the remainder functions lZ g ^i(z) collect 
terms that are independent of p^ to . 

Up to the order shown in the equations above, the functions F gg and entering the 
factorization formulae (13UI) and (1311) are closely related to corresponding quantities entering 
the resummed cross section for Higgs production at small transverse momentum [23J Ej- The 
only difference is that in that case L± = ln(x^/i 2 with 6 = 2e~ 7s . An exception is the 
two-loop coefficient d^ , which receives contributions from diagrams with two real emissions 
and is thus sensitive to the definition of the jet veto. 

The calculation of the kernel functions at one-loop order proceeds in close analogy to the 
calculation of the corresponding functions entering in the small-g^ resummation for Drell-Yan 
or Higgs production, which was performed in SCET in [23J HI]. The notation used above 
matches that adopted in these papers. The relevant Feynman graphs are shown in Figure [2j 
Since at this order only a single parton line is cut (corresponding to one-particle states X), 
the jet veto is trivially implemented as an upper cut on the transverse momentum of that 
parton. In addition to the collinear one-loop corrections shown in Figure [21 the corrections 
to the soft function in must be computed. At one-loop order this is trivial, since the 
relevant integrals are scaleless and therefore vanish in our regularization scheme. The results 
of our calculation agree with the above expressions, which were derived from the corresponding 
evolution equations. For the one-loop constant term in fl3^|) and the remainder functions in 
fl35l) . we obtain 



dr° = , 



o 



TZ g<r - q (z) = 2C F z . 



(37) 



The soft function S vanishes at one-loop order (like in the case of small-g^ resummation for 
Drell-Yan production [23]), but it is expected to be non-zero in higher orders, where the jet 
veto imposes non-trivial phase-space constraints. 

The two-loop coefficient d™ to in (J3IJ) cannot be determined from our one-loop analysis. 
It accounts for genuine two-loop effects and is therefore different from the corresponding 
coefficient d\ relevant for small-g^ resummation in Higgs production [231 EH] . The authors of 
[T7] have derived a partial expression for the /^-dependent piece of the jet-veto cross section 
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at two-loop order in perturbation theory, finding that 



a(p v f to ) 



00 0) 



= Q4C A a 

2— loop 



f(R) In 

„veto 
PT 



(38) 



where the dots represent -R-independent terms and other, possibly independent terms that 
are not enhanced by the large logarithm ln(m#/p™ to ). They observed that the function f(R) 
has a smooth behavior in the limit where R — > oo and adopted a normalization in which it 
vanishes like 1/R in this limit [45] . For very large R, the jet clustering algorithm combines 
all particles into a single jet, since d^ — > in (jlj). The jet veto then implies an upper bound 
on the vector sum of the transverse momenta of all particles, which is equivalent to a bound 
on the transverse momentum of the Higgs boson. Assuming that the formalism of [T7j is 
still valicfl in the limit of very large R, we can then match the -R-dependent two-loop terms 
obtained in this paper with the ones following from (l3Tj) and from the results obtained in our 
previous work [231 EI] for the transverse-momentum spectrum of the Higgs boson. In this way, 



we obtain 

or equivalently 



dr°=4 + 2C 3 (r£) 2 -8r£/(a), (39) 



d v 2 cto / 808 _ \ „ 224 



C A V 27 *°J 27 



— + 4C 3 \C A - — T F n f - 32f(R) , (40) 



where f{R) vanishes for R — > oo. For realistic values R < 1, this function can be very well 
approximated by the expansion [TJJ 

f(R) = - (1.0963 C A + 0. 1768 T F n f ) In R+ (0.6072 C A - 0.0308 T^n/) 

(41) 

- (0.5585 C A - 0.0221 T F n f ) R 2 + (0.0399 C A - 0.0004 T F n s ) R A + . . . . 

It would be desirable to check expression f)39p with an explicit calculation of the coefficient 
^veto j n gQET This is left for future work. We have pointed out in [23] that the coefficient d g 2 
entering in the small-gr resummation for Higgs-boson production is related to the correspond- 
ing coefficient d\ for Drell-Yan production by the Casimir-scaling relation d\/C F = d 2 /C A , and 
that in this ratio only maximally non-abelian color structures arise. Such a scaling relation no 
longer holds in the presence of a jet veto, because the function f(R) does not share this prop- 
erty. In particular, in the Drell-Yan case f(R) contains some "abelian" terms proportional to 
C F [17], which give rise to a two-loop C\ term in the exponent of the anomaly factor, indicating 
that this exponent is no longer constrained by the non-abelian exponentiation theorem. 

The explicit results collected in this section, along with the known NLO expressions in RG- 
improved perturbation theory for the hard matching coefficients C t in ([2]) and C$ in (1A4[) and 
[5]), provide the basis for the first NNLL resummation of the Higgs-boson production cross 



section with a jet veto. In the context of effective field-theory calculations it is conventional 
to count orders in RG-improved perturbation theory, while the literature on collider physics 
usually quotes the logarithmic accuracy of the resummation. A resummation with NNLL 
accuracy in the exponentiated cross section corresponds to a cross-section calculation at NLO 



1 See the note added in proof. 
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in RG-improved perturbation theory. At this order, we need the two-loop anomaly function 
in (I34p but only one-loop expressions for the kernel functions in (135]) . Before presenting the 
numerical results of such an analysis, it is interesting to compare our approach with that 
recently put forward in [17], where the first resummation of the jet- veto cross section at NLL 
order was accomplished. The results of these authors can be recovered from our expressions 
by using the lowest-order approximations 

F gg (p™ to , fi) = a s It L ± , I g ^(z, p™ to , fi) = 5(1 - z) 5 gi , (42) 

along with the leading-order RG-improved expressions for the Wilson coefficients C t (which 
equals 1 at this order) and Cs- Our results go significantly beyond this approximation, and 
they can in principle be extended to higher orders in a systematic way. We will comment 
below on the numerical differences between NLL and NNLL calculations for Higgs production 
at the LHC. 

While a two-loop calculation of the kernel functions I g ^i is beyond the scope of this work, 
we believe that such a calculation may still be practical using presently available tools. Com- 
bined with the known expressions for the hard matching coefficients Ct and Cs at NNLO 
in RG-improved perturbation theory, this would provide the most important ingredients for 
a calculation of the jet-veto cross section at N 3 LL order. (For a complete result, also the 
four-loop cusp anomalous dimension and the three-loop anomaly coefficient c?3 eto would be 
needed, but their numerical effect is expected to be rather small.) Sensitivity to the jet algo- 
rithm arises first in two-loop diagrams in which two intermediate lines are cut, corresponding 
to final states X containing two collinear, anti-collinear, or soft partons. In the phase-space 
integrals for these partons, we can change variables such that 

d 3 k { dyid^ihridhn . 1 „ , AO , 

— ; « = 1,2, (43) 



(2tt) 3 2^ 16tt 3 

where hr% = \k±i\, and without loss of generality we can set <p\ = and 02 = <P- If for the 
two emitted partons (yi — y 2 ) 2 + 2 < R, then these two partons are merged into a single 
jet, and the jet veto must be applied to the sum of their transverse momenta, i.e., we must 
impose that a/ k\ x + k\ 2 + 2A>ri cos <fi < p™ to . This happens inside a disk with radius R 
in the y-<f) plane. If the first condition is not satisfied, then the partons are classified as 
two individual jets and the jet veto must be applied to each of them, i.e., we must require 
that /ctij&T2 < Pr to - The integrals over transverse momenta can be evaluated after the jet 
clustering has been performed. 



4 Higgs production at the LHC 

We now perform a phenomenological analysis of the jet-veto cross section at the LHC with 
i/i = 8TeV. For concreteness, we consider the cross section integrated over rapidity, given 
in ( 13T|) . We set R = 0.4 and recall that to the order we are working the results are the 
same for all inclusive, kx-style jet algorithms (kr, anti-£>r, or Cambridge/ Aachen). The R 
dependence only enters via the coefficient d™ to in ( 139]) and the matching corrections, both of 
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which are NNLO terms which involve at most two real emissions. It follows that the clustering 
is the same for all version of the jet algorithm. All our numerical predictions are obtained 
using MSTW2008NNLO PDFs with a s (M z ) = 0.1171 [H]. While one could argue that it is 
more natural to evaluate the lower-order predictions with NLO PDF sets, we prefer to keep 
the PDFs fixed in order to asses the size of higher-order corrections. Furthermore, even the 
resummed predictions at NLL order contain the dominant higher-order corrections, so that 
using NNLO PDF sets seems more appropriate. 

We take into account finite quark-mass effects in the Born-level cross section and the 
leading electroweak corrections by replacing 



<T (p) ->■ <r (ji) k E w J ^ A(x q ) ; x q = — / - ie (44) 

q=t,b 



4m 2 



m H 



in (131]) . where 



x + x(l — x) arctan 2 



1 



(45) 



The value of this function is very close to 1 for the top quark, while A(xb) = 0(m 2 /m 2 H ) * s 
rather small. We use m t = 172.6 GeV and set the Higgs mass to ran = 125 GeV. For the 
6-quark mass, we use m&(m&) = 4.2 GeV, which translates to ra&(m.ff) = 2.83 GeV used when 
evaluating the cross section. When this is done, the absolute value of the sum |-4(xt)+A(x6)| = 
0.9997, indicating that quark-mass effects are tiny at Born level. We include electroweak effects 
as calculated in [16J H7] HH] and implemented in the code RGHlGGS jl9] for the total cross 
section. For m# = 125 GeV, this gives a factor kew = 1-0514. Overall, the change in the 
Born- level cross section obtained from (I44p amounts to an enhancement by a factor 1.0507. 

The first plot in Figure [3] shows our resummed predictions for the jet- veto cross section at 
NLL and NNLL order. These results are obtained by evaluating relation (13 ip using ( I34p and 
( 135]) as well as the NLO expressions for the hard matching coefficients C< and Cs obtained in 
RG-improved perturbation theory. For the default value of the factorization scale we adopt 
the canonical choice \x = p™ to , while the default hard matching scales are /it = m t in (T5]) and 
fi 2 h = —m 2 H in (IA4p . The bands in the figure are obtained by varying the three scales indepen- 
dently by a factor 2 about their default values and then adding the resulting uncertainties in 
quadrature. The main effect is due to the variation of the factorization scale \i. At NLO, the 
change in the cross section from the \Xh and fit variations is only about 1%. 

Next, we turn to the treatment of the hard function C 2 (m 2 ,n) \Cs{— /i)| 2 in the fac- 
torization formula (131 p . which arises as an overall factor both in the jet- veto cross section 
and the total cross section for Higgs production. As shown in [HI E], the Wilson coefficient 
Cs receives very large perturbative corrections, which arise in the analytic continuation of 
the scalar form factor Cs{Q 2 , Hh) to the time-like momentum transfer Q 2 = —m 2 H relevant 
for Higgs production. These corrections can be resummed by choosing the time-like value 
fi\ = —m 2 H for the hard matching scale and then using RG evolution to evolve back to space- 
like scale choices. The positive effect of this scale choice can be seen by comparing the first and 
second plots shown in Figure [3] The red bands in the left panel correspond to scale variations 
about our default choice /x 2 = —m 2 H , while the blue bands in the middle panel correspond to 



16 




Figure 3: Resummed predictions for the Higgs-boson production cross section with a jet veto 
at the LHC (with y/s = 8 TeV). The plot on the left shows our results at NLL and NNLL order 
obtained with the default choices of matching scales and scheme. The second plot illustrates 
the presence of the large corrections encountered when a space-like value of the hard matching 
scale Hh is used. The third plot shows the effect of factoring out the hard function at NNLO. 



variations about /x| = m 2 H . With such a choice, the resummed result suffers from the same 
large higher-order corrections that affect the fixed-order expansion of the total cross section. 

Since the Wilson coefficients Ct and Cs are known even at NNLO in RG-improved per- 
turbation theory, and since they appear as prefactors in the expression for the jet-veto cross 
section, it is not unreasonable to use the most accurate expressions for these coefficients that 
is available. Numerically, the difference between taking the NLO and NNLO hard functions 
is not very large. Including the NNLO corrections leads to a 2% increase of the cross section 
at low p^ to and a 1% increase at higher values. Using the NNLO approximation for the hard 
function will also be convenient below, when we match our resummed expression to fixed- 
order perturbation theory. From here on, we write the RG-improved hard function in the 
exponentiated form 

Ct(m 2 t ,/j,)\C s (-m 2 H ,fi)\ 2 = exp 2g 1 (m 2 , /i t , //) + 2Reg 2 (-m 2 H , fi h , //) , (46) 

where the functions g\ and g 2 will be expanded to NNLO in the exponent. The corresponding 
expressions are given in the appendix. The last plot in Figure |3] shows that this treatment 
improves the convergence behavior of the expansion for the resummed cross section, since the 
bands obtained at NLL and NNLL order now overlap. We stress, however, that while the 
scheme ambiguities reflected by the results in the figure are significant when the resummation 
is performed at NLL order, the differences between the various schemes become negligible 
once we perform the resummation at NNLL order. This demonstrates the importance of our 
analysis, which for the first time has enabled a resummation with this accuracy. 
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To obtain the best possible predictions for the jet-veto cross section, we match our re- 
summed expression (13TT) to the known NNLO fixed-order result for the cross section, which 
can be obtained by running the codes FeHiP [50] or HNNLO [ST] 152] . In our numerical 
analysis below we use the latter. The matching to fixed order allows us to also include terms 
that are suppressed by powers of p™ to /mH, which are not captured by our leading-power fac- 
torization formula. In addition to the power-suppressed terms, the matching also provides 
us with the two-loop constant, which is not predicted by a resummation at NNLL order. To 
ensure that the matching correction is not contaminated by the large perturbative corrections 
contained in the hard function, we will factor out the hard function (as discussed above) and 
define a reduced cross section cr(p^? to , /S) as 

a(p v T cto ) = C t 2 K,/i) \C s (-m 2 H ,fi)\ 2 x • (47) 

To remove the prefactor from the fixed-order result, we need to divide it by the two- loop fixed- 
order expansion of the hard function, which is given in flA8j) in the appendix. The matching 
is then performed for the reduced cross section according to 

onnll+nnlo = Snnll + (onnlo - ONNLL| cxpandcd to NNLO ) • (48) 

Our result for the resummed cross section (at NNLL order) matched to fixed-order pertur- 
bation theory (at NLO and NNLO) is shown in the first plot in Figure HI We observe that 
the matching corrections become more significant at larger values of the jet veto. This is 
expected, since the power corrections are controlled by the ratio p™ to /mH- Once the resum- 
mation is implemented in the way discussed above, the dominant matching corrections arise at 
NLO. The additional effect of performing the matching at two-loop order is to further reduce 
the residual scale dependence slightly. The second plot in the figure shows our results in a 
different scheme, in which in addition to the hard function also the anomaly term - the last 
factor in the first line of fl31|) - is factored out before the matching to fixed-order perturbation 
theory is performed. Comparing the two plots, we observe that the results are very similar, 
even though the scale variation at NNLL+NLO is somewhat larger in the second case. This is 
to be contrasted with the results shown in the last panel in Figure HI which are obtained if one 
neglects to factor out the hard function in the way shown in (1471) . In this case, the matching 
corrections suffer from a rather bad perturbative behavior, which becomes particularly obvious 
at small values of p™ to . 

The first plot in Figure [BJ shows our final results for the resummed and matched jet- 
veto cross section. For comparison, in the second panel we also give the result when the 
resummation is only performed at NLL order. After matching to NNLO, the difference between 
the two results is rather small but not entirely negligible. For example, at p^ to = 20 GeV the 
band obtained at NNLL+NNLO is approximately 4% higher than at NLL+NNLO. However, 
at NNLL order the matching corrections are much smaller. The excellent stability of the 
result is comforting, but is in part due to our improved matching scheme, in which the hard 
function is always taken at NNLO in RG-improved perturbation theory. Then the difference 
between the matched NLL and NNLL results is only due to higher-order terms from the 
anomaly. In Table [TJ we give predictions for the jet-veto cross section including estimates 
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Figure 4: Cross-section predictions obtained in the default matching scheme (left) and an 
alternative scheme (center), in which also the anomaly term is factored out before the matching 
to fixed-order perturbation theory is performed. The right plot illustrates the bad perturbative 
behavior encountered when the hard function is not factored out. 



of the perturbative uncertainty (shown by the bands in Figure [5]) as well as PDF and a s 
uncertainties. The latter are estimated by re-evaluating the cross section with a number of 
error PDFs, which were fitted to data for different a s values and encode the experimental 
uncertainties in the data sets entering the PDF fits [TT]. The uncertainty on the cross section 
is then computed from these results according to the prescription associated with a given PDF 
set. We find that the relative PDF and a s uncertainties are rather insensitive to higher-order 
effects. To excellent accuracy, they can be inferred from the Born-level cross section evaluated 
at the low scale \x = p^ to . It is thus a simple matter to compute the PDF uncertainty also for 
other PDF sets. 

We emphasize that in our approach a prediction is obtained directly for the Higgs-boson 
production cross section with a jet veto - the observable that hopefully will soon be measured 
at the LHC. In contrast, the authors of [17] focused on the jet-veto efficiency, defined as 
e(p v rf to ) = o"(p^ to )/cr tot . In order to compare to their results, we now briefly discuss the veto 
efficiency as well. For the total cross section, we use the central value of the result obtained 
in [H] and implemented in the code RGHlGGS At y/s = 8TeV and m H = 125 GeV, 

we obtain cr tot = (19-65io]iJi:L48) pb> where the second error gives the combined PDF and 
a s uncertainties at 90% CL. The resummed expressions for the jet-veto and the total cross 
sections are quite similar. They involve exactly the same hard function, and to NLL order 
the relative PDF and a s uncertainties are the same when evaluated at the same factorization 
scale. These uncertainties therefore cancel to a large extent when taking the ratio of cr(p™ to ) 
and crtot, and hence they can be neglected in the predictions for the jet- veto efficiency. The 
scale uncertainties in the two cross sections are correlated in such a way that the relative error 
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Figure 5: Resummed and matched predictions for the jet-veto cross section (first two plots) 
and efficiency (last plot) for Higgs-boson production at the LHC. 



on e(p™ to ) would be slightly smaller than that on <j(p™ to ). To be conservative, we normalize by 
the central value of the total cross section and retain the full scale uncertainty of the jet-veto 
cross section. The corresponding results are shown in the third plot in Figure as well as in 
the last column of Table [TJ To a good approximation, one can infer the jet- veto cross sections 
corresponding to other PDF sets by multiplying the veto efficiency with the central value for 
the resummed total cross section and add its PDF uncertainties. The corresponding values for 
cr tot can be obtained by running the RGHlGGS code with a different set of PDFs. For example, 
we find a tot = (20.87±g;H±i;t|) pb using NNPDF2.0 [S3], and a tot = (19.72±g;ff±?;^) pb using 
CTEQ6.6 PDFs P]. 

Our perturbative results for the jet-veto cross section have a high accuracy: the scale 
uncertainty from varying the factorization scale fi by a factor 2 around fi = p™ to is about 
12% for j9™ to = lOGeV and decreases to 4% at p™ to = 30GeV. This is almost as small as 
the scale uncertainty on the total cross section, which is 3%. The uncertainty from varying 
the hard matching scales fi t and fi^ is negligible, since we evaluate the hard function at 
NNLO in RG-improved perturbation theory. Indeed, by now also the three-loop correction to 
Cs(—m 2 H —ie, /i/i) is known [551 ESI ETJ , and it amounts to a mere —0.015% shift for /i^ = —m 2 H . 
(With = +m 2 H , the correction would be 4%, which confirms once again the advantage of 
using a time-like scale choice.) While the logarithmic accuracy of the resummation is one 
order higher for the total cross section, our result for the jet-veto cross section includes the 
most crucial part, the hard function, at the same accuracy, and so it should not come as a 
surprise that the scale uncertainties turn out to be comparable. 

Finally, let us relate our findings to the NLL+NNLO results presented in the recent paper 
|17j . An important difference to our work is that these authors do not compute the cross section 
directly, but instead present results for the jet-veto efficiency, arguing that its uncertainties 
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pvcto a(j0 v C to) [p b ] e ( p veto) 

in A Q7+0.59+0.19 n 9C -o+0.030 

1U <±.J< _o. 01-0. 22 u - zoo -0.001 

1 c; 7 QQ+0.62+0.34 n oQQ+0.032 

lO I •0-J_g.05_o.35 u - oyj -0.002 

9H in 1 Q+0.61+0.48 nt;ii;+0.(l31 

ZU 1U.10„ 08-0.49 u - oiJ -0.004 

25 11.88_ffiS 0.605_^ 

'in 1 OQ+0.48+0.71 n fi7 o+0.024 

OU lO.ZO„ 03 _o 70 U.U(O_ 02 

oo 19.6618^15* 1 



Table 1: Predictions for the NNLL+NNLO resummed and matched cross section and efficiency. 
The second error on the cross section gives the combined PDF and a s uncertainties at 90% CL. 

are largely uncorrelated with those of the total cross section. In that way they try to avoid 
the large perturbative corrections known to affect the gg — >■ H amplitude. Indeed, since 
the hard function is the same for the jet-veto and the total cross section, the corresponding 
large corrections formally cancel in the ratio of these two quantities. However, the extent 
to which this cancellation takes place depends on the precise way in which the perturbative 
expansion of the efficiency is performed. Equivalent matching schemes therefore lead to sizable 
differences in their results. Specifically, the authors use three matching schemes, related to 
different definitions of the fixed-order expansion of the efficiency e(p™ to ). In their scheme (a), 
the efficiency is obtained by expanding both the resummed jet- veto cross section and the fixed- 
order total cross section to a given order and then taking the ratio of the two. In scheme (c), 
the perturbative expansion is performed directly for the efficiency, while scheme (b) infers the 
efficiency from the normalized 1-jet rate. The fixed-order result for (a) and (c) is shown in the 
right two panels of Figure [TJ As we discussed in the introduction, the difference between these 
two schemes is so large because the fixed-order result for the jet- veto cross section involves two 
types of large corrections: those from the analytic continuation of the hard function, which 
also affect the total cross section, and those from Sudakov logarithms associated with the jet 
veto. Different schemes combine these in different ways, which then leads to quite different 
matching corrections. Even for p™ to = 20 GeV, the difference between schemes (b) and (c) is 
more than 20%, which is shocking given that it is formally an 0{a^ s ) effect. The origin of this 
problem is that the authors of [T7j resum the Sudakov logarithms associated with the jet veto, 
but they do not control the large corrections to the total cross section. 



5 Summary 

We have presented for the first time a factorization theorem for the Higgs-boson produc- 
tion cross section in the presence of a veto against jets with transverse momentum above 
a threshold p™ to . It forms the basis for the systematic resummation of Sudakov logarithms 
a" ln fc (m^/p™ to ) with k < 2n beyond the LL approximation. In analogy with the situa- 
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tion encountered in the analysis of the transverse-momentum spectrum of the Higgs boson, 
the final expressions ( 130]) and ( 13T|) are affected by a collinear anomaly, which introduces an 
extra power-like dependence on the hard scale m#. We have derived the all-order form of 
this anomaly and have determined its two-loop expression by combining information on the 
two-loop transverse-momentum spectrum of the Higgs boson with a recent result for certain 
jet-radius dependent two-loop terms obtained by Banfi, Salam and Zanderighi in [17] . Using 
results on the RG-improved hard matching functions available from previous studies of Higgs 
production in SCET, we have been able to perform the first resummation of the jet- veto cross 
section and efficiency at NNLL accuracy. These results have then been matched to the avail- 
able NNLO fixed-order result for the cross section, which accounts for power corrections in 
the ratio pff to /mH as well as the non-logarithmic two-loop contributions. 

In addition to Sudakov logarithms of the ratio of m#/p™ to , the cross section also suffers from 
large perturbative corrections to the hard function, which arise in the analytic continuation of 
the scalar gluon form factor to the time-like kinematics relevant for the gg — >■ H amplitude. 
The same corrections also affect the total cross section. They can be resummed in a systematic 
way using RG methods 0, [7J. To obtain the most accurate prediction possible, we have 
performed the matching of our NNLL result to fixed-order perturbation theory after the hard 
function has been factored out, as shown in (j47|) . Once this is done, the matching corrections 
are quite small, and we observe a very good perturbative stability of our results. 

Our phenomenological results presented in Section H] demonstrate that, with our improved 
formalism, we are in a position to predict the jet-veto cross section with high accuracy. We 
stress that we obtain the physically measurable cross section directly, and not by combining a 
prediction for the veto efficiency with a separate calculation of the total cross section. In fact, 
at least in one respect the theoretical analysis of the jet- veto cross section is simpler than that 
of the total cross section: it is evident that resummation is necessary and that non-logarithmic 
terms are strongly suppressed at small p™ to . The resummation of the total cross section, on 
the other hand, is based on an expansion around the partonic threshold, and it is not a priori 
clear how much of the cross section is generated by this region. Our final results for the cross 
section with a jet veto, shown in Figure [5] and collected in Tabled], exhibit scale uncertainties 
decreasing from 12% at p™ to =10 GeV to 4% at p™ to = 30 GeV. Additional uncertainties from 
PDFs and the value of a s are of a similar magnitude. The uncertainty we achieve is similar 
to the one for the total cross section. 

The results obtained in this paper form the basis for model- independent, high-precision 
calculations of the Higgs-boson production cross section in the presence of a jet veto. We be- 
lieve that, using available theoretical methods, it may be possible to extent the NNLL+NNLO 
analysis presented here to include the dominant logarithmic corrections arising at N 3 LL order. 
Even without such an improvement, however, the accuracy of our predictions is sufficiently 
high to allow for a reliable prediction of the cross section, which can hopefully soon be con- 
fronted with experimental data. 
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Note added in proof: After this paper was accepted for publication, reference [65] appeared, 
in which a NNLL resummation formula for the jet-veto cross section was derived that is 
equivalent to our result (}3~Tj) . However, in this work it was also shown that the result (1380 . 
which we have taken from [T7], is modified in the limit of large R. The additional contribution 
is such that the term proportional to £3 in (I39p gets cancelled. The correct relation then reads 
^vcto _ ^9 _ gY^f^R), and as a result in (HI)]) the coefficient of £3 in the first term becomes —28 
instead of 4. The numerical impact of this correction on our results will be studied elsewhere. 

In the recent paper |66j it was claimed that the clustering of soft and collinear particles 
into the same jet is unsuppressed for R = 0(1), and hence our factorization formula f[3~Tj) 
would break down. If correct, this observation would imply that neither our result nor that 
of [65] correctly reproduce the structure of the cross section at NNLL order. We have con- 
vinced ourselves that the soft-collinear mixing terms are absent if the multipole expansion is 
correctly implemented in SCET, or equivalently, if the double counting of overlapping mo- 
mentum regions is removed by performing appropriate subtractions. For the example of two 
independent emissions studied in [66] , we have explicitly rederived the result (4) of [63] from a 
SCET analysis based on our formula fl3TT) . We also note that the authors of [65] have checked 
their NNLL corrections against available fixed-order results up to 0(o? s ). 
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A Results for the hard matching coefficients 



The Wilson coefficient Ct in the effective Lagrangian (JI]) satisfies the RG equation [24] 



d 



C t (m 2 t , /J 2 ) = 7*(a a ) C t (m 2 , fj 2 ) , with 7 *(o; s ) 



d /3(a a 



da;. 



(Al) 



The anomalous dimension is related to the QCD /3-function /3(a s ) = da s /d In /i, because 
the two-gluon operator in the effective Lagrangian is the Yang-Mills Lagrangian. The exact 
solution to the evolution equation reads 



C t (m 2 ,y) 



C t (m 2 ,fi t ) , 



(A2) 



where [i t ~ m t is a hard matching scale at which the Wilson coefficient can be evaluated 
in fixed-order perturbation theory. The two-loop expression for C t (m 2 ,fi t ) was calculated in 
[58| 159] and is given in eq. (12) of [TJ. The three-loop expression is also known [601 l6Tj . 

The RG equation for the SCET matching coefficient of the gluonic current in ([6]) reads [7] 



d 



C s (-m 2 H - ie,n) 



cusp 



<iln/i 

and its general solution has the form 
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/4 



7T 



(A4) 
(A5) 



is the matching condition at an appropriately chosen hard scale /i^, at which the Wilson 
coefficient has a well-behaved perturbative series. We have demonstrated in [7J that this is 
not the case when one adopts the naive choice fih ~ n^H- On the other hand, a time-like scale 
choice such that fi\ ~ -m|j gives rise to a nicely convergent series. The fact that a s (fih) is 
complex in this case does not pose any difficulty for the numerical evaluation. The two-loop 
corrections to (1A5[) can be derived from the results of [H2] and have been given in eq. (17) 
of [7J. The three- loop corrections are also known [551 ESI [57J- In the solution (1A4f) we have 
introduced the definitions 1631 



o s (ii) 



T A a) 

da 



a a {v) 



(3(a) 



a s {v) 



da' 



a s (fj.) 



a T {v, At) 



T A (a) 

i cusp V ) 

da ■ 



a B (y) 



(3(a) 



(A6) 



and similarly for the function a 7 s . The perturbative expansions of these functions obtained at 
NNLO in RG-improved perturbation theory can be found in eqs. (93) and (94) in the appendix 
of [Mj. The required anomalous-dimension coefficients entering these expressions are given, 
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e.g., in eqs. (A2) of [7] and (95) of [M]. Up to three- loop order, the cusp anomalous dimension 
in the adjoint representation differs from that in the fundamental representation simply by an 
overall factor Ca/Cf- 

The functions g± and g 2 entering the exponent in (j46|) follow from the above results. We 
find 

/%s(/^))Awt) 

g 2 (-m 2 H , /j, h , n) = 2S(fi h ,n) - a r (jih,(*) m — r^" ~ V^/uAO + l n Cs(- m L A*h) • 

For the matching of our resummed result for the jet- veto cross section to the corresponding 
fixed-order expression, it is necessary to derive a fixed-order result for the reduced cross section 
a(p™ to , /i) defined in (14 7p . in which the hard function is factored out. To compute it to NNLO, 
we need to divide the fixed-order result by the two-loop expansion of the hard function. For 
N c = 3 colors and rif — 5 light quark flavors (the top quark is integrated out), we obtain 

CUmlri \Cs(-m 2 H ^)\ 2 = 1 + d ^ + c 2 (^) 2 + • • • , (A8) 

where 

ci = -6L 2 h + 22 + 7tt 2 

r4 46 „ / 698 ,\ r2 /1240 184tt 2 \ r 
c 2 = 18L* + -L\ + - 36tt 2 J L 2 h + f — — - 36C 3 J L h (Ag) 

274^ 10718 1679tt 2 998Cs 37tt 4 
~~ * + 27 + 6 3" + 2 

Here Lh = \n(m 2 H / fi 2 ) and L t = ln(m 2 //i 2 ) are the large logarithms. 
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